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ABSTRACT 

Photon-limited transit surveys in V band are in principle about 20 times more 
sensitive to planets of hxed size in the habitable zone around M stars than G 
stars. In I band the ratio is about 400. The advantages that the habitable zone 
lies closer and that the stars are smaller (together with the numerical superior¬ 
ity of M stars) more than compensate for the reduced signal due to the lower 
luminosity of the later-type stars. That is, M stars can yield reliable transit de¬ 
tections at much fainter apparent magnitudes than G stars. However, to achieve 
this greater sensitivity, the later-type stars must be monitored to these corre¬ 
spondingly fainter magnitudes, which can engender several practical problems. 

We show that with modest modihcations, the Kepler mission could extend its 
effective sensitivity from its current My = 6 to My = 9. This would not capture 
the whole M dwarf peak, but would roughly triple its sensitivity to Earth-like 
planets in the habitable zone. However, to take advantage of the huge bump in 
the sensitivity function at My = 12 would require major changes. 

Subject headings: planets 


1. Introduction 

The one conhrmed transiting planet discovered to date, HD209458b (Gharbonneau et al. 
2000), lies only 10 R© from its host G star, and therefore well inside the so-called “habitable 
zone”, where water could exist in its liquid state. All ongoing transit surveys are primarily 
sensitive to such “hot Jupiters” because their large diameters give rise to relatively strong 
(~ 1%) signals, while their proximity to their host increases the probability of transits 
occurring (Howell et al. 2000; Brown & Gharbonneau 1999; Mallen-Ornelas 2002; Udalski et 
al. 2002; Burke et al. 2002; Street et al. 2000, 2002). 
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While gas giants are not expected to themselves be habitable regardless of their loca¬ 
tion, they could have water-laden moons like Europa, Ganymede, and Callisto, which could 
support life if the planet lay in the habitable zone. Future ground-based surveys could be 
sensitive to such habitable-zone giant planets. Moreover, the planned Kepler mission^ has 
as its primary goal the detection of planets down to Earth size in or near the habitable zone. 
The Kepler targets would be even more direct analogs of the Earth. 

Here we investigate the sensitivity of transit surveys to planets in the habitable zone 
as a function of stellar type. We derive two remarkable conclusions. First, transit surveys 
are most sensitive to habitable planets orbiting middle M stars {My ~ 12), despite the fact 
that they are most sensitive overall to planets orbiting early G stars (My ~ 4). Second, 
to achieve this sensitivity to dim stars, the survey’s magnitude limit must be an extremely 
strong function of spectral type. Specihcally, it must be about 7 magnitudes fainter for 
middle M stars than G stars. Since such variable requirements can signihcantly affect survey 
design, they deserve careful analysis and consideration. For example, we show that unless 
the Kepler mission extends its magnitude limit for M stars well beyond its current (type 
independent) limit oi V = 14, it will lose most of its potential sensitivity to habitable 
planets. 


2. Sensitivity to Habitable Planets 


Pepper, Gould, & DePoy (2002) showed that if every star of luminosity L and radius R 
has a planet of radius r in a circular orbit of semi-major axis a, then a photon-noise limited 
transit survey will detect 
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planets. Here n{L, R) is the number density of stars of the specihed type, G is the angular 
area of the survey, and dg is the distance out to which an equatorial transit can just be reliably 
detected for the hducial parameters Lg, i?g, rg, and og. The numerical factor rj ~ 0.719 
arises because the volume probed by the survey is smaller by a factor for non-equatorial 
transits, where y is the ratio of the transit chord to the stellar diameter. For the most part we 
will be comparing the detectability of planets of the same size around different stars. Hence, 
except when we are evaluating the normalization of equation (1) for a specihc example, we 
will ignore the last factor. 


^http://www.kepler.arc.nasa.gov/ 
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The equilibrium temperature of a planet (and thus whether or not water can exist in 
liquid phase on the planet surface) is the same if a oc where the bolometric luminosity 
Lboi is to be distinguished from L, the luminosity in the band of observation. Hence, the 
relative sensitivity of a transit survey to habitable planets is 

N oc (2) 

Here we assume that all planets have the same albedo and neglect atmospheric effects. If we 
compare, for example, G stars {My = 5) with middle M stars {My = 12), the ratios of the 
various factors are Ni2/N^ ~ 6 X 630-3/2 ^ 47/2 ^ gQ5/4 10. Note that if we were comparing 

detectability at the same semi-major axis rather than the same habitability, the last factor 
would not have entered, and the ratio would have been 0.05. Hence, while the sensitivity to 
planets in general is completely dominated by G stars, the sensitivity to habitable planets is 
completely dominated by M stars. That is, the lower luminosities (and so smaller semi-major 
axes) combined with the smaller radii and greater numerical density of M stars more than 
compensate for the reduced photon counts. 

In Figure 1, we show the sensitivity to habitable “Earths” and to Earths all at the 
same semi-major axis (taken to be 1 AU). That is, the histograms show the total number 
of planets N that will be detected as a function of My assuming that each star in the held 
has one Earth-size planet in the habitable zone or, respectively, one such planet at 1 AU. 
(The two histograms cross at My = 5 because the habitable zone is then at 1 AU). The 
absolute normalization of this plot is set according to the characteristics of the Aep/er mission 
(7.8 X lO^e-hr-^ at U = 12, G = 105deg2, Ay = 0.3, mission-total S/N =8 required for 
detection), but the form of the histogram would be the same for any photon-limited survey. 
The normalization for any other planet size should be multiplied by a factor (r/r^)^, and 
the normalization for any other hxed semi-major axis should be multiplied by (a/AU)-3/2 
(see eq. [1]). The hgure is constructed assuming that detection is in V band. The effect of 
substitnting other bands is approximately to change the slope of the histogram. For example, 
since the slope of the main seqnence is dMy/d{V — I) = 3.37 (Reid 1991), snbstitntion of I 
band would lead to an increase of slope dA\ogN/dMy = (3/2) x 0.4/3.37 = 0.178. That is, 
middle M stars wonld gain relative to G stars by an additional factor of = 18. 

To compute these histograms, we follow the procedure of Pepper, Gould, & DePoy 
(2002). An important feature of the color-magnitude diagram is that while the main sequence 
is fairly narrow for My > 6 , it broadens for brighter stars (due to faster stellar evolution), 
so that a star of a given My can have a large range of colors. Thus, for the well-dehned 
lower main seqnence. My > 6 , we consider the luminosity function (Bessell & Stringfellow 
1993; Zheng et al. 2001) in 1 mag bins, and evalnate the stellar radius at the center of each 
bin using the linear color-magnitnde relation of Reid (1991), the color/snrface-brightness 
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relation of van Belle (1999), and the VIK color-color relations from Bessell & Brett (1988). 

On the other hand, for the upper main sequence My < 6 we evaluate the histograms 
directly using the Hipparcos catalog (ESA 1997). For example, the luminosity function for 
My = 4 is computed by summing X^j[(4/3)7rZlf]“^ over all stars within the Hipparcos com¬ 
pleteness limit, V < 7.3, having 3.5 < My < 4.5, and lying within 50 pc. The distance Di 
is the minimum of 50 pc and the distance at which the star would have V = 7.3. Then, the 
constant-semi-major-axis histogram is computed by summing (and appropriately normal¬ 
izing) L^^^i?~'^^^[(4/3)7rZlf]“^, while the habitable-zone histogram is found by summing 

The stellar radii are determined from Hipparcos/Tycho 
{Bt,Vt) photometry and the color/surface-brightness relation of Gould & Morgan (2003), 
ultimately derived from van Belle (1999). We evaluate Lboi using bolometric corrections as a 
function of H — / color derived from Binney & Merrifield (1998) and Bessell & Brett (1988) 
at the bright end and Reid & Gilmore (1984) at the faint end. 


3. Magnitude Limits for Habitable Planets 

To achieve the sensitivities calculated in § 2, one must analyze the light curves of 
all the stars being probed. This statement would appear so obvious as not to be worth 
mentioning. However, as we now show, dim stars are being probed at substantially fainter 
apparent magnitudes than are their more luminous cousins. Hence, to avoid losing most of 
the sensitivity of a transit experiment one must set different magnitude limits for stars of 
different My. 

Let miiin be the apparent magnitude at which the survey achieves the minimum accept¬ 
able signal-to-noise ratio (S/N) for a planet of radius r and semi-major axis a circling an 
My = 5 star. Gonsider an identical planet circling an My = 12 star in the same orbit and at 
the same apparent magnitude. Since the star’s radius is a factor 4 smaller, the S/N will be 
a factor 4^ x 4 “^A = 43/2 |;ijxies larger. The first factor comes from the fact that the planet 
occults a larger fraction of the stellar surface, the second from the reduced total time spent 
in transit (“duty cycle”) due to the smaller radius. Now, the bolometric luminosity of the 
dimmer star is down by a factor 80, so to keep in the habitable zone, the planet must be 
moved closer by a factor 80 ^A. This increases the transit duty cycle by the same factor and 
so increases the S/N by 80^/^. Gombining these two effects with the usual dependence of 
S/N on flux implies S/N oc Hence, to maintain the same S/N, one must 

increase the magnitude limit by. 


Amiim = 0.5AMboi — 7.5Alogi?. 


(3) 
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For the example just given, Amum = 6.9, i.e., a factor 570 in flux. Note that Amum does not 
depend on the passband of observation. The dashed curve in Figure 1 shows the difference 
in the maximum distance modulus relative to the Sun at which a habitable planet (of given 
radius) can be detected. Specihcally, A(m — M)o = (Vjim — My) — (Vjim,© — My^o). For V 
band observations, this offset is roughly constant over 10 magnitudes. That is, the limit¬ 
ing magnitude required to achieve the sensitivities shown by the bold histogram increases 
approximately in lock step with My. This behavior can have major consequences for the 
design of transit experiments. 


4. Discussion 

Photon-limited observations are not routinely achievable over very large dynamic ranges. 
A number of effects can intervene at fainter magnitudes. First, once the stars fall below the 
sky, the number of detected systems no longer falls as (see eq. [1]) but as L^. That is, 
once the sky is reached, detections would fall off by a factor ~ 4 per magnitude relative to 
the bold histogram in Figure 1. Second, if the exposure times are set so as not to saturate 
the brightest target stars (say 2 mag brighter than the mag limit appropriate for My = 5), 
then the flux levels will be a factor ~ 4000 times lower at the mag limit appropriate for 
My = 12. For many observing setups, this would render the flux levels comparable to the 
read noise. Third, a held that is very uncrowded at one magnitude may well be extremely 
crowded 7 magnitudes fainter, and this may interfere with doing photon-limited photometry. 
Hence, it is unlikely that all of the potential shown in Figure 1 can be achieved in any given 
practical experiment. Nevertheless, this potential is so great that it is worth thinking about 
how to achieve as much as possible. 

We illustrate the type of trade-offs involved by considering the Kepler mission, the only 
transit experiment proposed to date whose primary target is planets in the habitable zone. 
The Kepler design calls for recording only stars V < 14, presumably in order to minimize 
data transmission. The mag limit required to detect Earth-like planets around Sun-like stars 
at S/N=8 is Wm,© = 13.6. Hence, from the dashed curve in Figure 1, the required mag limit 
at My = 6 is Viira = Wm,© “ My^Q + My + A(m — M)o = 14.4, i.e., already somewhat fainter 
than the V = 14 limit adopted by Kepler. The volume probed, and hence the detections, fall 
relative to the bold histogram by a factor (10*^'^)^ ~ 4 for each magnitude beyond My = 6. 
See Figure 2. Hence, assuming that every star has one planet in the habitable zone, only a 
total of 20 will be detected^. 


^Actually, we have not taken account of stellar variability in this analysis, nor of binaries. Each of these 
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Keeping the photometry information for stars to V = 17, for example, should increase 
the number of habitable-planet detections by almost a factor 3, while keeping stars with 
K < 20 would increase this number by a factor 9. See Figure 2. Of course, there are many 
obstacles to keeping information on stars this faint. 

For example, there would be an enormous number of giant stars contaminating such 
a deep sample. In fact, however, distinguishing the later-type dwarfs from the much more 
numerous giants of similar colors is straightforward. As shown by Gould & Morgan (2003), 
they can easily be identihed on a reduced proper motion (RPM) diagram constructed using 
data from USNO-A (Monet 1996, 1998) and 2MASS (Skrutskie et al. 1997). With USNO-B 
(Monet et al. 2002) it should be possible to extend coverage to fainter magnitudes and to 
achieve higher precision as well. The total number of such stars My < 12 is only ~ 3 x 10^, 
far fewer than the ~ 10^ giants in the held that could be eliminated using the same RPM 
diagram (Gould & Morgan 2003). 

A more difficult problem is crowding. The Kepler point spread function (PSF) is de¬ 
liberately defocused to 10", meaning that the crowding limit is about lOarcmin”^. At the 
adopted Aep/er sight line, (/,6) = (69.6,5.7), this limit is reached at Rusno ~ 18-2 (as de¬ 
termined by a query of the USNO-A catalog). If Kepler were redirected to (/,6) = (70,15), 
the density of stars down to Rusno = 18.2 would be 3.5arcmin“^, approximately 3 times 
lower. Hence, recovery of faint stars would be much easier. Presumably, Kepler has chosen 
to look right in the Galactic plane because of the higher overall star density. However, since 
the stars useful for transits mostly he within 500 pc, a held at b = 15° would he only 130 
pc above the plane, where the density of such target stars is only slightly lower than at the 
plane. Hence, the crowding problem could be substantially ameliorated at small cost. 

The last problem is sky. At high ecliptic latitude, the sky in space is U ~ 23.3 arcsec“^, 
or U ~ 17 in a 10" PSF. Hence, Kepler completeness appears to be fundamentally limited 
to stars My < 9. As shown in Figure 2, detections from stars at My > 10 would be highly 
suppressed. Only contracting the PSF (which has been deliberately defocused to improve 
the photometry) could overcome this difficulty. (Note that at U ~ 17, crowding would not 
be a serious problem even for the current Galactic-plane line of sight.) 

While the above remarks apply specihcally to Kepler, any transit experiment attempting 
to detect habitable planets would face similar constraints and trade-offs. 


might plausibly reduce detections by 25%, the former by making transits unrecognizable, the later by making 
habitable planets dynamically unstable and by increasing the total light in the aperture. Hence, the true 
detection rate for this optimistic scenario should be reduced by roughly 50%. 



Finally, we note that the results reported here are a strong function of planet size. As 
remarked in § 2, the absolute normalizations in Figure 1 scale as r®. By an argument similar 
to the one given in § 3, Amum = 10A log r. For example, for r = 0.63 r®, the normalizations 
in Figure 1 would be reduced by a factor 16, while the required Vnm at each My would be 
brighter by 2 mag. Hence, in Kepler^s current conhguration, it would retain full sensitivity to 
My ~ 8 stars, and so would be able to detect ~ 3 planets (see Fig. 2), while if its coverage 
were extended 3 mag to Viim ~ 17, it could probe the M-dwarf peak, and so detect ~ 8 
planets. 


We thank Scott Gaudi for his detailed comments on the manuscript. This work was 
supported by grant AST 02-01266 from the NSF. 
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Fig. 1.— Total number of planets detected in the habitable zone {bold histogram) and in 
1 AU orbits {solid histogram) as functions of absolute magnitude My (assuming every star 
has one such planet). The absolute normalizations have been set to Earth-size planets and 
to the characteristics of the Kepler mission, but the form is completely general assuming 
observations in V band. In / band, the slope would be tilted upward toward faint stars. 
Extinction has been taken into account, but stellar variability and binarity have not. The 
calculation uses different methods for early {My < 6) and late {My > 6) stars. The solid 
and dashed histograms compare the two methods in the one bin of overlap. The dashed curve 
shows the difference (relative to Sun-type stars) in the limiting distance modulus to which 
stars of different My must be monitored to achieve the sensitivity shown by the “habitable 
zone” histogram. Since this curve is basically flat, dimmer stars must be observed to the 
same distance and hence to much fainter apparent magnitudes. 
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Fig. 2.— Total number of planets detected normalized to characteristics of Kepler and 
under various assumptions. Bold solid histogram: each star has one Earth-size planet in the 
habitable zone and all stars are monitored regardless of magntiude (same as Fig. 1). Bold 
dashed histogram: same planet distribution, but now only stars V < 14 are monitored. Thin 
dashed histogram: same planet distribution, V < 17 stars are monitored. Dotted histogram: 
same planet distribution, V < 20 stars are monitored. Thin solid histogram: planets have 
radius r = 0.63 r® and only V < 14 stars are monitored. 













